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Abstract. It is shown that inhomogeneous Szekeres and Stephani universes exist 
corresponding to non-dissipative binary mixtures of perfect fluids in local thermal 
equilibrium. This result contradicts a recent statement by Zarate and Quevedo (2004 
Class. Quantum Grav. 21 197, Preprint gr-qc/0310087), which affirms that the only 
Szekeres and Stephani universes compatible with these fluids are the homogeneous 
Friedmann-Robertson- Walker models. Thus, contrarily to their conclusion, their 
thermodynamic scheme do not gives new indications of incompatibility between 
thermodynamics and relativity. Two of the points that have generated this error are 
commented. 
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1. Thermodynamic scheme for a mixture of two perfect fluids 

In a recent paper, Zarate and Quevedo extend the standard thermodynamic scheme, 
of local thermal equilibrium for a simple perfect fluid, to the case of a mixture of two 
perfect fluids. 

They deflne their thermodynamic scheme by the following relations: 

(a) conservation of the perfect fluid energy-momentum tensor: 

V-T = 0, T = {p + p)u^u + pg, (1) 

where p is the total energy density, p the mixture pressure and u the 4-velocity of the 
total matter flow. 

(b) conservation of the total current density: 

V • (nu) = , (2) 
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where n is the total particle number density, 
(c) Gibbs thermodynamic relation: 

T ds = d{p/n) + pd{l/n) — ^dc , (3) 

where s is the entropy per particle, T the temperature, c the fractional concentration of 
one of the two perfect fluid components and is the mixture chemical potential, i.e. the 
difference between the chemical potentials of the components of fractional concentration 
c and 1 — c respectively. 

Under these relations, and as Zarate and Quevedo emphasize, the entropy 
production no longer vanishes but is a consequence of the change in the fractional 
concentrations of the components. 

In addition, this thermodynamic scheme for a mixture of two perfect fluids recovers 
the standard scheme for a simple fluid when c = constant. In this sense their scheme is 
more general than the standard one. 

Nevertheless, its physical applicability risks to be much more restrictive than that 
of a simple fluid: the physical applicability of the non-dissipative (relativistic) evolution 
hypothesis clearly diminishes generically when in addition c ^ constant, i.e. the fluid 
is submitted to endothermic or exothermic reactions. A more realistic model would be 
consisted in a mixture of two perfect fluids generating a semi-perfect fluid, that is to 
say, a dissipative Pascalian one, with eventually a heat flux proportional to the gradient 
of the fractional concentration. 

Anyway, as no complete set of constraints is known on the space of formal 
constitutive relations in continuous thermodynamics (i.e. the more or less general 
inequalities on the thermodynamic variables and on some of their partial derivatives 
are insufficient to separate physical from unphysical equations of state), it makes sense, 
as Zarate and Quevedo emphasize, to analyze the compatibility of their thermodynamic 
scheme with Einstein equations. 

2. Compatibility of Zarate and Quevedo's thermodynamic scheme with 
Einstein equations 

In the article P] that motivates this comment, the authors quote a work by us P] 
and they assert: "Coll and Ferrando have shown that an exact solution [to Einstein 
equations] admits a thermodynamic scheme provided the integrability conditions of 
Gibbs equation are satisfied" . 

It is here pertinent to observe that, of course, this fact is the starting point of the 
quoted paper, but that our main goal there is i) to obtain these integrability conditions 
explicitly and in terms of the sole hydrodynamic variables {u,p,p), and ii) to interpret 
them physically. Our corresponding result is [2j (see also |3j): i) a divergence-free perfect 
fluid energy tensor evolves in local thermal equilibrium if, and only if, the space-time 
function ^ = p/p depends only on the variables p and p: dC, /\ d p A dp = 0; ii) then ^ 
is a state variable, C, = ^{p,p), and represents the square of the velocity of the sound. 
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As, on one hand, Einstein equations are biunivocally related to the sole 
hydrodynamic variables (m, p,p) and, on the other hand, the condition d^AdpAdp = 
is not a consequence of them, it follows that: there exist perfect fluid solutions to Einstein 
equations that do not admit a standard thermodynamic scheme, i.e. that cannot he 
interpreted as evolutions, in local thermal equilibrium, of a single fluid. 

What is the corresponding result for Zarate and Quevedo's thermodynamic scheme 
for a non-dissipative mixture of two fluids? 

As there exists only one scalar constraint for the compatibility between Einstein 
equations and the relatively restrictive standard thermodynamic scheme (namely, the 
one expressing that the quotient p/p is a function of state: p/p = C,{p,p)), it seems 
that the corresponding result for Zarate and Quevedo's scheme would be the absence 
of constraint, because of its less restrictive character. But this appears as a startling 
evaluation facing Zarate and Quevedo's main statement in [1 . Let us see it in detail. 

Suppose given, in a domain of the space-time, a solution to Einstein equations for 
a perfect fluid, {u,p,p) (in fact, a solution to the divergence-free equations dH), the 
argument that follows being also valid for test fluids). Then, u, p and p are known 
quantities in that domain and, in particular, the conservation equation for the total 
current density Q becomes a linear and homogeneous scalar equation in the scalar 
unknown n. Pick in it a particular solution n; with these four known elements {u, p, p, n) 
we can evaluate the one-form uj defined by 



Then, according to Gibbs equation (jHl) for the mixture, one has to explore the existence 
of the four thermodynamic scalars (s, c, T, p) such that 



But this existence is exactly what Pfaff decomposition theorem locally guarantees in 
four dimensions for any one-form uj , and a fortiori for the restricted one given by (jl)) ; 
we have thus: 

Proposition 1 Any perfect fluid solution to Einstein equations is compatible with a 
Zarate and Quevedo's binary thermodynamic scheme. 

How many Zarate and Quevedo's schemes may be associate to any given perfect 
fluid solution to Einstein equations? 

From the data {u,p,p) obtained as a solution to equation (jT]), the solution to 
the density conservation (j2)) is obviously determined up to an w-invariant function / , 
/ = u°'daf = , so that n = fuo , Uo being a particular solution. And for every 
such n , equation (jH) determines the particular one-form oj . Then, equation (0) defines 
generically the four thermodynamic variables s , functions of the three 

ones p , p , n ; the functions s and c result thus involved by one first order differential 
equation, namely 




(4) 



uo = T ds + pdc 



(5) 




(6) 
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which always admits solutions in one of the unknowns for every arbitrary choice of the 
other, such solutions depending (for example via an initial problem) on an arbitrary 
function of two of the three variables. Then the variables T and /x are univocally given 
by 



where 



r s 



h = n{s'/p - s'pc'p) . (8) 



Proposition 2 The different Zdrate and Quevedo's binary thermodynamic schemes that 
any given perfect fluid solution {u,p,p) to Einstein equations admits, are obtained by 
the free choice of a u-invariant function (determination of the total particle number 
density n), of a function of three variables (say, the entropy per particle, s{p,p,n)) 
and of a function of two variables ( say, the fractional concentration at a given value k 
of n, c{p,p,k)). 

The main statement in pjj asserts that among the Szekeres and Stephani 
families of perfect fluid cosmological models, the only ones that admit a [Zarate and 
Quevedo's] binary thermodynamic scheme are the Friedmann-Robertson-Walker models. 
Proposition 1 shows that this statement is wrong. Proposition 2 gives an indication of 
the 'distance' between this statement and the correct one. 

As already said, no complete constraints are known on the space of formal 
constitutive relations insuring the physical character of a model; nevertheless, the 
richness of the choice of Zarate and Quevedo's schemes stated in Proposition 2 locally 
guarantees the usual thermodynamic inequalities (such as T > or < c < 1). 



3. Some Zarate and Quevedo's thermodynamic schemes for Szekeres and 
Stephani universes 

The Szekeres and Stephani universes admitting a standard thermodynamic scheme, i.e. 
that of a one-component perfect fluid, have been considered by different authors. Bona 
and Coll |lj have shown that the Stephani universes admitting a thermodynamic scheme 
are space-times with a 3-dimensional group of isometrics acting on 2-dimensional 
orbits. This result has been recovered in [S] where the authors have also studied 
the thermodynamic Szekeres-Szafron models, and they have shown that a family of 
thermodynamic solutions of class II without symmetries exists. A different approach 
has been used in [Hj in order to study the Szekeres-Szafron space-times of class II: 
the solutions which represent a perfect fluid in local thermal equilibrium have been 
obtained and the associated thermodynamic schemes explicitly given. More specific 
thermodynamic analysis have been considered in other works. Thus, Sussman |B] has 
presented a family of Stephani universes which admit, up to a good approximation, an 
interpretation as classical mono-atomic ideal gases or as matter-radiation mixtures. On 
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the other hand, we have obtained all the Stephani universes which represent a generic 
ideal gas in local thermal equilibrium [7j. 

It is worth pointing out that, although in some cases the standard thermodynamic 
scheme imposes symmetries on the metrics, these papers show that inhomogeneous 
Szekeres and Stephani universes exist and are known, that admit a standard 
thermodynamic scheme. 

On the other hand. Propositions 1 and 2 show that all perfect fluid space-times 
admit Zarate and Quevedo's thermodynamic schemes, but they give no explicit solutions 
to them. The difliculties to flnd such explicit solutions associated with a general perfect 
fluid lie only in the solution of the two differential equations ((21) and (jHl) because, as we 
have seen, then the entropy per particle s (or equivalently the fractional concentration 
c ) may be chosen arbitrarily, and the temperature T and the mixture chemical potential 
fi are explicitly given by ((Zj) and 

Nevertheless, it is not difficult to flnd explicit Zarate and Quevedo's thermodynamic 
schemes for any perfect fluid admitting a standard thermodynamic one. Remember that 
a standard thermodynamic scheme is a perfect fluid solution {u,p,p) to equation (jT)) 
that admits a solution n to equation and for which there exist functions s and T 
of, say, the variables p and n , verifying 

Tds = d{p/n) +pd{l/n) . (9) 

Let {u, p,p,n,s,T) be such a standard thermodynamic scheme. Choose a space- 
time function c = c(a;") satisfying < c(x") < 1 , but otherwise arbitrary, let $(c) be 
an arbitrary real function of c , and deflne a new function s of p , n and c by 

s = s{p,n,c) = s{p,n) + <!>{c) , (10) 
Then, on account of Q one has: 

Tds = d{p/n) +pd{l/n) +T^\c)dc (11) 
and, consequently, calling 

/i = -T$'(c) , (12) 

the set of variables {u, p, p, n, s, T, c, p) deflnes a Zarate and Quevedo's thermodynamic 
scheme. 

Consider then the Szekeres and Stephani universes that admit a standard 
thermodynamic scheme and that are explicitly given in the above mentioned references 
[3-7]. On every one of them, every choice of a pair of functions c(a;") and $(c) directly 
deflnes, by means of equations (|TO|l and (fT2|l . a Zarate and Quevedo's thermodynamic 
scheme, i.e. an explicit counterexample of the main statement of 

4. Remarks 

Perhaps it is worthwhile to indicate a pair of incorrect arguments in Zarate and 
Quevedo's work P that could explain their erroneous conclusion. 
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In dealing with Szekeres universes, they select, by integration of a space-like 
equation, the sole solution /i = {p + p)/n , neglecting an arbitrary function of time (see 
their equation (27) in IJ). They neglect it under the argument that /i being a function 
of three thermodynamic variables, and the expression /x = {p + p)/n already containing 
them, one cannot add to fi an arbitrary function of time without increasing the number 
of independent variables. This argument is incorrect because it mix up considerations 
on space-time variables and thermodynamic ones without taking care of their specific 
connection: the pressure being a function of the sole time in Szekeres space-times, to add 
a function of time is nothing but to add a function of the pressure p , an addition which 
is perfectly admissible in their situation but that invalidates their particular expression 
of Gibbs equation (their equation (28) in [T]) and their consequences. 

In dealing with Stephani universes, they impose to them two equations of state, 
for fj, and s (their equation (40) in P), obtained in the Szekeres case under i) the 
abusive restriction above mentioned and ii) the hypothesis of spatially homogeneous 
pressure, p,i = , identically satisfied in Szekeres universes but generically inadmissible 
in Stephani ones. 
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